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ABSTRACT 

Starting from the observation that kilohertz Quasi Period Oscillations (kHz QPOs) occur in a very- 
narrow range of X-ray luminosities in neutron star low mass X-ray binaries, we try to link the kHz 
QPO observability to variations of the neutron star magnetospheric radius, in response to changing mass 
inflow rate. At low luminosities, the drop off of kHz QPOs activity may be explained by the onset of 
the centrifugal barrier, when the magnetospheric radius reaches the corotation radius. At the opposite 
side, at higher luminosities, the magnetospheric radius may reach the neutron star and the vanishing of 
the magnetosphere may led to the stopping of the kHz QPOs activity. If we apply these constraints, the 
magnetic fields of atoll (B ~ 0.3 - 1 x 10 8 G for Aql X-l) and Z (B ~ 1 - 8 x 10 s G for Cyg X-2) sources 
can be derived. These limits naturally apply in the framework of beat frequency models but can also 
work in the case of general relativistic models. 

Subject headings: stars: neutron — stars: individual (Aql X-l, 4U 1820-30, Cyg X-2) — X-ray: stars 
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1. INTRODUCTION 

Color-color diagrams (CDs) have been widely used to 
gain insight on low mass X-ray binaries (LMXRBs) host- 
ing a neutron star (e.g. Hasinger & van der Klis 1989). 
These sources are usually classified as atoll or Z sources 
depending on the path that they trace in the CDs as their 
luminosity vary. Atoll sources are characterised by an up- 
ward curved branch (lower and upper banana state) and 
one or more harder island states (cf. Fig. 1). Z sources 
produce a Z-shaped path in the CD with three branches, 
from hard to soft: horizontal branch (HB), normal branch 
(NB) and flaring branch (FB). Type I X-ray burst proper- 
ties correlate well with position in the CD (van der Klis et 
al. 1990). This led to argue that position on the CD along 
the atoll was a good indicator of M (but X-ray flux was 
not). For Z sources (a subset of which show sporadic X- 
ray bursts too), optical and UV emission as well as quasi 
periodic oscillations (QPOs) models indicate that M in- 
creases as HB^NB^FB (cf. van der Klis 1995). 

RossiXTE observations have shed new light on these 
sources, revealing QPOs in the kHz frequency range (for 
a review see van der Klis 2000). In atoll sources, kHz 
QPOs occur in the island and lower banana part of the 
CD (Mendez 1998 and references therein). In the case of 
Z sources kHz QPOs are present all the way down the HB 
and disappear in the lower part of the NB (e.g. Wijnands 
& van der Klis 1998b and references therein). KHz QPOs 
are therefore detected in a narrow range of fluxes and in 
similar positions in the CDs. Whereas the disappearance 
of the kHz QPO activity at the low luminosity end may be 
related in some atoll sources to a lack of sensitivity, this is 
certainly not the case at the bright end. 

In this letter we try to link the source behavior in the 
CD to the physics of the accretion process. We mainly 
concentrate on atoll sources, since for Z sources data are 
more sparse and they still lack an indication of the spin 
period from e.g. type I X-ray bursts. 



2. ACCRETION REGIMES 

Accretion phenomena onto a fast spinning, weakly mag- 
netic neutron star three radii which do not vary with the 
mass inflow rate can be identified: 

• the neutron star radius, R, which depends on the 
neutron star mass and equation of state. We scale 
here R = 10 R§ km for a neutron star of mass 
M = 2M 2 M Q ; 

• the marginally stable orbit radius, r ms . In the case 
of a rapidly rotating neutron star a small correc- 
tion to the Schwarzschild value has to be included. 
As a first order in the specific angular momentum 
j [j = 2ircI/(G M 2 P), with c the speed of light, 
/ and P the neutron star moment of inertia and 
spin period and G the gravitational constant] , r ms ~ 
(6 - 2.31 j)GM/c 2 (Kluzniak 1998). In the case of 
'soft' equations of state this radius is larger than R, 
whereas for stiff equation of states it is smaller and 
does not play any role. In the case of a spin fre- 
quency of 1.8 ms (as in the case of Aql X-l; Zhang 
et al. 1998a), and for a 2 Mq neutron star with 



I = 2 x 10 45 g cm- 2 , r n 



16 M 2 km; 



• the corotation radius, r cor , at which a test parti- 
cle in Keplerian orbit corotates with the neutron 



star, 
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In the case of a spin 
1.8 Pi. 8ms ms) we have 



As the mass inflow rate changes the limiting radius at 
which the neutron star magnetic field pressure balances 
the ram pressure of the mass inflow towards the neutron 
star (i.e. the magnetospheric radius, r m ) changes too. In 
the case of spherical accretion with a dipolar magnetic field 

4 \l/7 

fm ~ I „ „t r »v, (where M is the mass inflow rate onto 

1 G M M z / v 



2 



□ 



□ 



□ 



□ 



D 

red 



p 



-E3- 



□ 



□ 



0.6 



o 

o 
o 

K 



0.4 





1.6 



1.1 



2.2 



soft Color- 
ed a 



Fig. 1. — Color-color diagram of Aquila X— 1. The soft and hard colors are defined as the ratio of count rates in the bands 3.5 — 6.0 kcV 
and 2.0 — 3.5 keV, and 9.7 — 16.0 keV and 6.0 — 9.7 keV, respectively. Points in the banana branch are marked with closed circles, in the 
lowest island state with filled squares, in the middle island with open circles and in the extreme island with open squares. Typical error bars 
in the banana and the island states are shown. Black and grey symbols indicate segments with and without kHz QPOs, respectively. Taken 
from Reig et al. (2000). 



the neutron star magnetosphere and fi = B R 3 is the mag- 
netic moment). A parameter £ is usually considered to ac- 
count for the case of disk accretion, with £ ~ 0.5 — 1 (Ghosh 
& Lamb 1979, 1992). Wang (1995) pointed out that £ 
is in the range 0.87-0.95. In the following, we assume 
£ = 1. For an accretion luminosity of L = GM M/R = 
10 37 L37 erg s -1 , we derive a magnetospheric radius of 
r m = 20 Bl /7 Rf /7 M\ n L~i' 7 £ km (where B = B$ 10 8 
G is the magnetic field). For higher luminosities the radi- 
ation pressure starts playing a role (region A in the stan- 
dard disk model) and the magnetospheric radius changes 

to r£ = 33 Bf 39 L 3 - 2/13 M 2 1/39 R^ 13 £ km (cf. Ghosh & 
Lamb 1992; here we neglect the weak dependence on the 
viscosity parameter). 

3. APPLICATION TO ATOLL SOURCES 

As a working example we can consider the case of 
the transient atoll source Aql X-l. This is the pro- 
totype of soft X-ray transient sources (Campana et al. 
1998a) and, together with the accreting X-ray binary pul- 
sar SAX J1808.4-3658 (Wijnands & van der Klis 1998a; 
Chakrabarty & Morgan 1998), is the only LMXRB which 
showed signs of centrifugal inhibition of accretion (Cam- 
pana et al. 1998b; Zhang, Yu & Zhang 1998; Gilfanov et 
al. 1998). During a type I burst of Aql X-l observed with 
RossiXTE, a coherent periodicity was discovered at 549 
Hz (Zhang et al. 1998a). This frequency may correspond 



to the spin frequency of the neutron star or, like the case 
of 4U 1636-53, to twice its value (Miller 1999). In the 
following we adopt a spin period of 1.8 ms. 

In Fig. 1 we show the CDs of Aql X-l based on 
RossiXTE. Data fall into three main distinct groups, which 
can be readily identified with the extreme island, island 
and banana states based on the hard color. By adopting 
an absorbed blackbody plus power law model, Reig et al. 
(2000) derived the 2-10 keV X-ray luminosity in the dif- 
ferent states of the CD diagram (for a distance of 2.5 kpc). 
This varied from 2.0 x 10 34 erg s _1 in the extreme island 
state (not shown in Fig. 1) to 3.2 x 10 35 erg s _1 in the 
island state. In the banana state, the luminosity increased 
from left to right, from 2.4 x 10 36 erg s _1 to' 3.9 x 10 36 
erg s _1 (Reig et al. 2000). Assuming the spectrum at a 
comparable luminosity observed with BcppoSAX (Cam- 
pana et al. 1998b), we derive a bolometric correction of 
~ 4 in the 0.1-40 keV energy band. We remark, however, 
that it has long been recognised that the observed flux 
does not provide a confident measure of the accretion lu- 
minosity. Different contributions may come from e.g. the 
neutron star rotational energy which is extracted by the 
disk (Priedhorsky 1986). 

Despite the large variations in X-ray luminosity, kHz 
QPOs are observed in selected regions of the CDs. In 
the transient sources Aql X-l and 4U 1608-52 (Mendez et 
al. 1999) there exist RossiXTE observations in which kHz 
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QPO activity was not detected. The corresponding lim- 
iting RossiXTE PCA rates for the two sources are ~ 600 
and ~ 400 c s" 1 , respectively (2-16 keV; Mendez 1998). 
Moreover, the persistent atoll source 4U 0641+091 showed 
kHz QPOs in the range 250-1800 c s^ 1 (2-60 keV; van 
Straaten et al. 2000). We conclude that, at least in the 
case of Aql X-l, the disappearance of the kHz QPO activ- 
ity at the faint end is likely not due to poor statistics. For 
what concerns the maximum flux at which kHz QPOs have 
been detected, no problems of statistics exists (~ 1500 and 
- 3000 c s" 1 for Aql X-l and 4U 1608-52, respectively, cf. 
Mendez 1998). The stopping occurs for all the atoll sources 
in the lower banana state. 

The existence of selected regions in the CD (even if 
slightly different from source to source) in which kHz 
QPOs are detected is noteworthy. We now try to iden- 
tify the onset and the disappearance of kHz QPO activity 
with changes in the magnetospheric radius. 

(i) As the accretion rate increases, the onset of the kHz 
activity can be related to the exit of the source from the 
propeller regime, when most of the matter is inhibited to 
accrete onto the neutron star surface by centrifugal forces 
(i.e. r m ~ r cor ) and therefore kHz QPOs can be hardly de- 
tected. The luminosity jump expected for the transition 
from the propeller to the accretion regime in the case of a 
ms spinning neutron star is just a factor of a few (Corbet 
1996; Campana & Stella 2000). This flux reduction is in- 
dependent on the mechanism invoked to explain the kHz 
QPOs, like beat-frequency models (BFMs; Strohmeyer et 
al. 1996; Miller, Lamb & Psaltis 1998) or general relativis- 
tic models (CRMs, Stella & Vietri 1998, 1999). 

The minimum luminosity for which kHz QPOs have 
been detected during the fading end of an Aql X-l out- 
burst is 1.1 x 10 36 erg s" 1 (2-10 keV; Zhang et al. 1998a). 
This luminosity is a factor of <~ 2 higher than the lumi- 
nosity corresponding at the steep increase in the hard- 
ness ratio and associated to the centrifugal barrier onset 
as interpreted by Zhang, Yu & Zhang (1998) and Cam- 
pana et al. (1998b). Considering the bolometric correc- 
tion given above we assume a 0.1-40 keV luminosity of 
L m in = 4.2 x 10 36 erg s _1 . This luminosity represents an 
upper limit to the accretion luminosity (which defines the 
magnetospheric radius), but at the onset of the centrifu- 
gal barrier external contribution such as the luminosity 
extracted from the neutron star rotational energy should 
weaken considerably. 

Equating the magnetospheric to the corotation radius, 
we derive an upper limit to the magnetic field value of: 

B < 1 x 10 8 L\H R- 5 ' 2 Ml" P 7 / 8 6 ms r 7/4 G. (1) 

In the case of Miller's et al. BFM the relevant radius for 
kHz QPOs is the sonic radius, r s , which lies inside the mag- 
netospheric radius. Since, r s < r m < r cor , the constraint 
above applies as well. 

(ii) The disappearance of the kHz QPOs at high lu- 
minosities may be related to the disappearance of the 
magnetosphere, when r m becomes comparable to the neu- 
tron star radius, R. A similar line of reasoning has been 
applied by Cui et al. (1998) interpreting the transition 
from island and banana and the disappearance of the kHz 
QPO, as the disk becoming disangaged from the magnetic 
field. In particular, by comparing r m with r ms Cui et al. 



(1998) were able to limit the neutron star magnetic field 
to - 0.6-0.9 x 10 8 G. 

The disappearance of the kHz QPOs can be interpreted 
in the framework of BFMs as the drop out of the mag- 
netosphere without which the beating mechanism cannot 
work. In the Miller's et al. BFM the relevant radius is the 
sonic point radius which, by definition, lies always inside 
the magnetospheric radius. On the other hand, the disap- 
pearance of the magnetosphere may leave no space for an 
orbiting blob to acquire the necessary eccentricity to gen- 
erate kHz QPOs as predicted by CRMs (Stella & Vietri 
1999) or if the magnetosphere itself is needed to generate 
the small eccentricity in a resonant interaction (Vietri & 
Stella 1998). In general both the sonic point model and 
the GRMs do not require the presence of a magnetosphere 
to work. However, we expect that the neutron stars in 
LMXRBs do have a magnetic field if they are the pro- 
genitors of millisecond radio pulsars. This link has been 
recently confirmed by the detection of coherent pulsations 
in the transient LMXRB SAX J1808-3658 (Wijnands & 
van der Klis 1998a). 

In the case of Aql X-l, kHz QPOs stop to be detected 
at a luminosity higher than 3.6 x 10 36 erg s _1 (2-10 keV; 
Zhang et al. 1998a). Given the small range of luminos- 
ity we adopt the same bolometric correction, resulting in 
0.1-40 keV luminosity of L max = 1.4 x 10 37 erg s _1 . In 
this case, we are going to derive a lower limit on the mag- 
netic field and no problems arc posed by taking a higher 
luminosity than the true accretion one. We can obtain a 
conservative estimate of the magnetic field from r m £ R: 

B > 3 x 10 7 L W R-"* M- 1 " r 7/4 G. (2) 

Different parameters can play a role in changing the exact 
value of these constraints as well as the uncertainties in 
the theory underlying the location of the magnetospheric 
radius. Despite these uncertainties, the two limits are well 
consistent and constrain the magnetic field of Aql X-l to 
B ~ 0.3 - 1 x 10 8 G. 

3.1. 4U 1820-30 

A further case is represented by the atoll source 4U 
1820-30 for which indication of a saturation in the kHz 
QPO frequencies versus the flux intensity has been ob- 
served (Zhang et al. 1998b; Kaaret et al. 1999 and refer- 
ences therein). 4U 1820-30 follows the general behaviour 
of atoll sources described above (Zhang et al. 1998b). The 
saturation starts near the lower banana onset (Zhang et 
al. 1998b). KHz QPOs appear at a 0.3-40 keV luminosity 
of 5 x 10 37 erg s _1 (for a distance of 8 kpc) and disappear 
at 8 x 10 37 erg s _1 , i.e. spanning a vary small range of lu- 
minosities. We adopt these as bolometric luminosities and 
derive, as described above, B £ 9 x 10 8 G (for a neutron 
star spin period of 3.6 ms, as inferred from the kHz QPO 
separation) and B £ 0.8 x 10 8 G, respectively. A further 
constraint comes from the interpretation of the saturation 
of the kHz QPO frequencies as due to the magnetospheric 
radius reaching the last stable orbit. If this is the case, the 
upper kHz QPO frequency is related to the Kcplerian fre- 
quency at the disk inner edge. In the case of GRMs, the 
ceasing of the kHz QPOs derives from the lack of space 
between the point at which the blob leaves the accretion 
disk and the neutron star. 
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A saturation value of 1065 Hz has been measured in 4U 
1820-30, leading to an estimate of the neutron star mass 
of ~ 2.1 Mq. Alternatively, this implies a radius of the 
marginally stable orbit of r ms = 17.4 km. The saturation 
occurs at a luminosity of 6.6 x 10 37 erg s^ 1 (0.3-40 keV), 
which provides a direct estimate of the magnetic field value 
of B ~ 2 x 10 8 G (close to the lower limit derived above 
as expected due to the fact that r ms ~ R). In the case 
of Miller's et al. BFM the estimate of the magnetic field 
gives only a lower limit, since the sonic radius is smaller 
than the magnetospheric radius. 

The saturation in the kHz QPO frequencies has been 
observed only in 4U 1820-30. A possible explanation re- 
lies on the large neutron star mass inferred for this ob- 
ject, which in turn implies a smaller radius than for a 
~ 1.6 - 1.8 M Q neutron star by a factor of 10-20%. This 
might be enough to prevent the observation of the kHz 
QPO saturation in lighter neutron stars. 

4. APPLICATION TO Z SOURCES 

In Z sources kHz QPOs are usually observed in the up- 
per part of the NB down to the end of the HB (e.g. Cyg 
X-2, Wijnands et al. 1998; GX 17+2, Wijnands et al. 
1997). Flux determinations are very rare. Only in the 
case of GX 349+2 (Zhang, Strohmeyer & Swank 1997) 
the onset of kHz QPO occurs at a 2-10 keV luminosity of 
~ 1.2 x 10 38 erg s _1 (for a distance of 8.5 kpc). This is 
consistent with the idea that the lower apex of the Z path 
occurs when the source attains the Eddington luminosity. 
Despite the variations observed in the CDs, Z sources vary 
very little in luminosity, with typical RossiXTE PCA rates 
in excess of 1000 c s _1 (Wijnands & van der Klis 1998b) 
such that problems with a too poor statistics do not exist. 

We consider the case of Cyg X-2, which showed kHz 
QPO down to apex NB/FB. The 0.1-20 keV luminosity ex- 
trapolated by EXOSAT observations at the NB/FB apex 
is Imax ~ 2 x 10 38 erg s" 1 (at 8 kpc; Chiappetti et al. 
1990). The luminosity is high and the inner disk regions 
are likely in the radiation pressure regime such that the 
magnetospheric radius should be adopted. Using this 
luminosity to infer the lower bound on the magnetic field 



we obtain B ^ 1 x 10 8 G. This is a conservative lower limit 
since at these luminosities the coupling between disk and 
neutron star should be high. KHz QPOs have been ob- 
served at the end of the HB, for which the minimum lumi- 
nosity observed with EXOSAT is just a factor of 2 smaller. 
From this we derive B ;$ 8 x 10 8 G, assuming a neutron 
star spin period of 2 ms, still waiting for the discovery of a 
coherent periodicity from bursting Z sources. Higher fields 
may be obtained by assuming faster spin periods. 

5. DISCUSSION 

We interpret the onset and the end of kHz QPOs in 
terms of the variable extent of the neutron star magneto- 
sphere. The turn-on of kHz QPOs at low fluxes is inter- 
preted as the starting of accretion onto the neutron star 
surface, overcoming the effects of the centrifugal barrier 
(here we do not claim that this barrier is working perfectly, 
but if it can stop the large fraction of the infalling matter, 
the flux reduction will prevent to detect kHz QPOs). The 
disappearance of the kHz QPOs corresponds instead to the 
'disruption' of the magnetosphere on the neutron star (see 
also Cui et al. 1998). Given these assumptions, we esti- 
mate in the case of the atoll source Aql X-l a magnetic field 
of B ~ 0.3 — 1 x 10 s G. Similar results can be obtained for 
other atoll sources. Interpreting the saturation of the kHz 
QPO frequency for increasing luminosities in 4U 1820-30, 
leads to a magnetic field estimate of ~ 2 x 10 8 . In the case 
of Z sources we derive a rough estimate based on luminos- 
ity obtained with EXOSAT. In the case of Cyg X-2, we 
obtain a very conservative lower limit of B ~ 1 — 8 x 10 8 
G. This is consistent with the idea that Z sources have 
higher magnetic field than atoll sources. 

This scenario allows for clear predictions on the lumi- 
nosities at which the kHz QPOs set in and out, which can 
be verified with RossiXTE monitoring of low mass X-ray 
binaries. 

I thank P. Rcig for permission of publishing Fig. 1. I ac- 
knowledge useful discussions with T. Belloni, M. Mendez 
& L. Stella. 
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